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ABSTRACT 


The launch vehicle is idealized as a one-dimensional structure consisting of a main 
beam to which several flexible branches are attached. The normal modes of both the un- 
restrained main beam and the cantilevered branches are determined. An energy ap- 
proach is then employed in which the displacement of the vehicle is expressed as the su- 
perposition of a finite number of free-free normal modes of the main beam and canti- 
lever normal modes of the branches. Lagrange’s equations are then used to derive the 
equations of motion in matrix form, and an iterative method of solution is included for 
completeness. 



METHOD FOR DETERMINING NORMAL MODES AND FREQUENCIES OF 
A LAUNCH VEHICLE UTILIZING ITS COMPONENT NORMAL MODES 

by James D. McAleese 
Lewis Research Center 

SUMMARY 

An analytical method is presented for determining the natural frequencies and lat- 
eral mode shapes of an idealized launch vehicle, which is envisioned as a one- 
dimensional structure consisting of a main beam to which several flexible branches are 
attached. Engines and sloshing masses are also included. The vehicle may be free or 
constrained by translational and rotational constraining springs. 

An energy approach is employed in which the displacement shape of the system is 
expressed as the superposition of a finite number of free-free normal modes of the main 
beam and cantilever normal modes of the branches. The choice of these modes as as- 
sumed displacement functions provides a simple, straightforward method of incorporating 
branches, engines, and sloshing masses in the analysis; many terms of the energy ex- 
pressions are eliminated and those terms that remain are easily evaluated. Lagrange T s 
equations are then utilized to develop the equations of motion in matrix form, and an 
iterative method of solution is included for completeness. 


INTRODUCTION 

In some instances, the launch vehicle can be appropriately idealized as a nonuniform 
Timoshenko beam for purposes of studying its lateral dynamic characteristics. In other 
cases, this model may not realistically represent the launch vehicle system. For in- 
stance, when large structural components are cantilevered within other portions of the 
vehicle, it becomes necessary to include them as flexible branches. An example of this 
is a spacecraft enclosed within a protective nose fairing, or an upper-stage engine sup- 
port structure extending into the interstage adapter. The inclusion of these branches 
becomes especially significant when clearance problems are anticipated or when the dy- 
namic loads within the branch itself must be evaluated. It may also be desirable to in- 



elude the inertial effect of the engines, which are generally represented as constrained 
rigid bodies, and the effect of sloshing propellants. 

The lateral normal modes of the launch vehicle system have been determined by a 
variety of methods. For example, in references 1 and 2, the main launch vehicle is 
analyzed as a nonuniform Timoshenko beam by the Myklestad method and the finite- 
element force method, respectively. In references 3 and 4, the vehicle is analyzed as 
a nonuniform Timoshenko branch beam system including engine effects by the Stodola 
method and the finite- element stiffness method, respectively. 

In the present analysis, the effects of flexible branches, swiveling engines, slosh- 
ing masses, and constraining springs can be incorporated into the system in a straight- 
forward manner. The normal modes of the unrestrained main beam and the normal 
modes of the cantilevered branches must first be determined. The displacement of 
the vehicle is then expressed as the superposition of rigid body displacement modes, 
normal modes of the free-free main beam, and cantilever normal modes of the branches. 
The energy expression is then formulated, and Lagrange’s equations are applied to de- 
rive the differential equations of motion. The principle of utilizing the uncoupled nor- 
mal modes of the system is not new; in reference 5, Scanlan proposes its application to 
airplanes, and Hurty (ref. 6), presents a general procedure for incorporating component 
modes in the finite- element displacement approach. The new feature in this analysis is 
the choice of this particular set of assumed modes to describe the normal modes of a 
launch vehicle system. This choice is logical when the nature of the connection of the 
branches in an actual launch vehicle is considered. 

The present approach is attractive because the normal modes of the main beam or 
of the branches can be developed by the most suitable method and then easily incorpo- 
rated in the system in terms of their modal characteristics. For example, a payload 
consisting of trusses and concentrated masses might be analyzed by the finite-element 
stiffness method and then coupled to the system by use of its cantilever normal modes, 
natural frequencies, and generalized masses. Also, describing a component in terms 
of its dynamic modes allows dynamic test data to be incorporated in the analysis. 


SYMBOLS 

A cross-sectional area of launch vehicle 

A effective shear area 

s 

(A g G) shear stiffness of i** 1 beam 

distance from attachment point of i beam to its center of gravity (i > 0) 
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q 0R 

q 0T 
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distance from center of gravity of system to center of gravity of main beam 

Young’s modulus 

bending stiffness of 1 beam 

total number of generalized coordinates used to describe motion of system 
number of elastic generalized coordinates used to describe motion of i beam 
shear modulus 

area moment of inertia of cross section 

mass moment of inertia of engine on r beam about center of gravity of engine 

i-T_ 

mass moment of inertia of i tn branch about its attachment point (i > 0) 
mass moment of inertia of main beam qbout its center of gravity 

stiffness of rotational spring at gimbal point of engine on i tn beam 

11. 

stiffness of r tn sloshing spring 

1 L. 

stiffness of j n translational constraining spring 

11. 

stiffness of j n rotational constraining spring 

total number of branch beams (does not include main beam) 

ir 

displacement of r m sloshing mass relative to equilibrium position 
displacement of r n sloshing mass relative to tank walls 

iu iu 

modal moment at clamped end of i tn branch for n cn uncoupled mode 
(i = 1> 2, . . . , k; n = 1, 2, . . . , f.) 

IV 

mass of engine on i in beam 
mass of i^ beam 

xi_ 

mass of r l sloshing mass 

IV IV 

modal shear at clamped end of i tn branch for n in uncoupled mode 
(i — 2 ? • • • , k, n - 1, 2, , f 

iu XI- 

generalized coordinate describing participation of n in uncoupled mode of i m 
beam in free vibration of system (i = 0, 1, . . . , k; n = 1, 2, . . . , fj) 

generalized coordinate describing rigid body rotation 

generalized coordinate describing rigid body translation 

kinetic energy of launch vehicle 

potential energy of launch vehicle 
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x longitudinal coordinate measured relative to center of gravity of system 

x gi coordinate of gimbal point of engine on i^ beam 
coordinate of attachment point of beam (i> 0) 
x r coordinate of r sloshing mass 

11. 

X/p j coordinate of connection point of j translational constraining spring 

J.V 

x 0j coordinate of connection point of j in rotational constraining spring 
y.(x) lateral displacement function of i th beam relative to equilibrium position 

11. 

angular rotation of engine on i beam relative to bending slope at gimbal point 

ii. 

j3j angular rotation of engine on i beam relative to equilibrium position 

e i distance from gimbal point of engine on i beam to center of gravity of engine 

IV 

e i (x) bending slope function of i beam relative to equilibrium position 
X^ displacement of center of gravity of swiveling engine on i^ beam 

M in generalized mass of n th uncoupled mode of i th beam (i = 0, 1, . . . , k; 

n = 1, • • • , f t ) 

ii. 

ju_ generalized mass of s L mode of vehicle 

p mass density 

IV 

pA i mass per unit length of i in beam 

IV 

plj mass moment of inertia per unit length of i cn beam 

(/?i n (x) displacement function of n th uncoupled mode of i th beam (i = 0, 1, . . . , k; 
n = 1, 2, . . . , fj) 

^ 0 r( x ) rigid body rotation displacement function 
<j9qt( x ) rigid body translation displacement function 

i^ in (x) bending slope function of n th uncoupled mode of i th beam (i = 0, 1, . . . , k; 
n = 1, 2, . . . , f.) 

^ 0 r( x ) rigid body rotation slope function 
i//q T (x) rigid body translation slope function 

w in natural frequency of n th uncoupled mode of i th beam, rad/sec 
(i = 0, 1, . . . , k; n = 1, . . . , fp 

co natural frequency of s th mode of vehicle 

s 

Matrices: 

[A] stiffness matrix 
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[B] inertia matrix 

[I] unit matrix 

{ q } column matrix of generalized coordinates 

Subscripts: 

i beams of system (i = 0, 1, 2, . . . , k) (main beam i = 0, branch beam i > 0) 

n, p uncoupled modes of beam (n = 1, 2, . . . , f^) 

Superscripts: 

’ differentiation with respect to time 

f 

differentiation with respect to x 

OUTLINE OF THEORETICAL DERIVATION 

The displacement of the launch vehicle is expressed as the superposition of the 
rigid body translation and rotation of the system, and a finite number of free-free modes 
of the main beam and cantilever modes of the branches. The amplitude of each of these 
assumed modes is governed by a generalized coordinate. The engines, which can gen- 
erally be considered as rigid in comparison with the supporting structure, are repre- 
sented as rigid bodies pinned at the gimbal point and constrained by rotational springs, 
as shown in figure 1. If the dynamic effects of sloshing propellants are to be included, 
the usual practice (ref. 7) is to introduce a mechanical analogy (spring mass or pendu- 
lum) into the system. In this analysis, sloshing is represented as a spring mass. For 
the case of a constrained launch vehicle, translational and rotational constraining 
springs may be attached to the main beam. In all, the total number of generalized co- 
ordinates is equal to the number of assumed modes plus the number of engines and 
sloshing masses. 

The kinetic and potential energies of the system T and U, respectively, are then 
expressed in terms of the assumed deflections and generalized coordinates. In the de- 
velopment of the energy expressions, both the main beam and the branches are treated 
as Timoshenko beams, but this is merely for ease of notation and does not affect the 
final equations. The fact that the assumed elastic modes are normal modes of specific 
components, and, hence, satisfy the orthogonality relations, eliminates many terms 
from the energy expressions. The remaining terms consist of physical and modal char- 
acteristics of the components; that is, mass, length, generalized mass, generalized 
forces, and natural frequencies. Applying Lagrange's equations yields the required dif- 
ferential equations of motion, the number of which is equal to the number of generalized 
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Actual Idealized 



Figure 1. - Launch vehicle configuration. 



coordinates. The equations in matrix form appear as 


[B] { 4 } + [A] { q } = {0} 

If harmonic motion is assumed, this equation is reduced to the eigenvalue equation 

m - [Ar l t B iV<i) - f o} 

which is then solved for the natural frequencies and associated mode shapes by an iter- 
ative method. 


Coordinate System 

As shown in figure 2, the launch vehicle configuration is described relative to a 
fixed x, y- coordinate system. The longitudinal coordinate x is measured relative to the 
center of gravity of the system. 

The lateral displacement is assumed to be representable as the superposition of 
rigid body modes and uncoupled normal modes of the components. Following this ap- 
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proach, the deflection of the main beam is expressed as the superposition of a rigid 
body translation mode <Pq t = 1; a rigid body rotation mode = x; and the free-free 
normal modes of the main beam cp qj, cp Qg, . . . , ; that is, 


y 0 (x,t) = <^ 0x (x)q 0T (t) + <P 0 R (x)q OR^ + ^01 (x)q 01 (t) + • • • + W 


= ‘POn^OnW d) 

n=T, R, 1 


where the coefficients 


q 0T’ q 0R’ q or * * * » q 0fQ 


are generalized coordinates that de- 


scribe the participation of the rigid body modes and free-free modes in the natural mode 
of vibration of the vehicle. Similarly, the bending slope of the main beam d q is 


6 q ( x > t) - , /' , Q'j'( x ) c lQ'j'(t) + ^OR^^OR^ + 1 ^01^ q 01^ + • • • + 


= Z , W x >W t > ( 2 ) 

n=T, R, 1 


where 

slope associated with translation mode <Pqj, 

slope associated with rotation mode <pQj^ 

bending slope functions associated with free-free modes of 
main beam 

The slope <^QT’ associated with the translation mode, is zero but it is carried along in 
the analysis for simplicity in notation. 

In like manner, the deflection of i tn branch y^ is expressed as the superposition of 
two rigid body modes associated with the translation and rotation of the branch as a 
whole, and the cantilever modes of the branch cp .^, cp^, • • • > <P^ ; that is, 

f i 

yj(x, t) = y 0 (x., t) + {x - Xjje^x., t) + (i = 1, 2, . . . , k) (3) 

n=l 


*0T =0 
^OR = 1 

^01’ ^02> • * * ’ ^0f 0 
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where Yq(Xj, t) is the displacement of the main beam at the attachment point x = of 
the branch, 0 q(x^, t) is the bending slope of the main beam at the attachment point, and 
q il’ q i2’ • ' ’ > q if are g enera li z ed coordinates describing the participation of the 

cantilever modes in the natural mode of vibration of the vehicle. When equations (1) 
and (2) are substituted into equation (3), the deflection is 

f 0 

y i< x > *> - [« , 0n< x l> + < x - *i > <'On< x l>] W‘> 

n=T, R, 1 


+ 2 ^in^ x ) q in (t ) (i = 1, 2, . . . , k) (4) 
n=l 

in 

Similiarly, the bending slope of the i in branch 6^ is 

f i 

6i(x,t) = 0 o (x.,t) + ^^ in (x)q in (t) (i = 1, 2, . . . , k) (5) 

n=l 

where • • • , i/'jf are the bending slope functions associated with the canti- 

lever modes of the branch. Substituting equation (2) into equation (5) gives the bending 
slope of the i Ln branch 

f 0 f i 

e i( x > t) = S ’W x i ) W t > + 2 ^in^in^ (i = 1, 2, . . . , k) (6) 

n= 1 n= 1 

After the deflection and bending slope of the main beam and branches are defined, the 
displacements of the swiveling engines and the sloshing masses can be defined. In the 
case of the engine, the displacement is expressed as the superposition of the motion of 
the clamped engine (the slope of the engine axis is equal to the bending slope at the gim- 
bal point) and the vibratory motion relative to the clamped position. The displacement 
of the center of gravity of an engine is 


A i (t) = y^Xg., t) + ejS^Xg., t) + e^t) (i = 0, 1, . . . , k) 


( 7 ) 


,th 


where y.(x t) is the displacement of the i beam at the gimbal point (x = x .), 0.(x t) 
l ci ei i ei 
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is the bending slope at the gimbal point, is the distance from the gimbal point to the 
center of gravity of the engine, and is the generalized coordinate describing the an- 
gular rotation of the engine relative to the bending slope at the gimbal point. 

The displacement of the center of gravity of a rigid engine can be expressed in 
terms of the generalized coordinates by substituting either equations (1) and (2) or equa- 
tions (4) and (6) into equation (7), depending on whether the engine is attached to the 
main beam or a branch. Therefore, for an engine on the main beam, 

f 0 

x 0^ = X) IW x eO> +e 0^0n( x e0 ) ] q 0nW + € 0 O! 0 (t) < 8 ) 

n=T, R, 1 


while for an engine on a branch, 
f 0 

\(‘>= £ [«’On( Jt l) + < x ei + £ i- x i) , W x i>] < W t > 

n=T, R, 1 


f i 


+ XKn( x ei> + e i^in( x ei)] q in^) + e i a i (t) (i = 1, 2, . . . , k) (9) 

n=l 

Following this same pattern, the angular rotation of the engine is 


/3.(t) = e^Xg., t) + o^t) (i = 0, 1, . . . , k) 


( 10 ) 


Substituting equation (2) into equation (10) gives the angular rotation of an engine on the 
main beam as 


- X ^On^eO^On^ + 
n=T, R, 1 


( 11 ) 


For the case of an engine on a branch, substitute equation (6) into equation (10) 


f i 


^i(l) = X) *0rS x i>%n® + X ^in^eiKn^ + (i = 1, 2, . . . , k) 

n=T, R, 1 n=l 


( 12 ) 
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The displacement of the sloshing mass L r can be written in terms of yQ(x r , t), 
the displacement of the main beam at the attachment point, and -£? r (t), the generalized 
coordinate describing the displacement of the sloshing mass relative to the attachment 
point, 


Ij r ( t ) = y 0 ( x r> *) + ( 13 ) 

Substituting equation (1) into equation (13) gives the displacement of the r 1 * 1 sloshing 
mass 


L r (t) = 2 ^0n( x r) q 0n W + ^ 

n=T, R, 1 


(14) 


Kinetic Energy 

The use of Lagrange’s equation requires the development of the kinetic energy ex- 
pression in terms of the generalized velocities. For the idealized system of beams, 
masses, and engines, the kinetic energy is 


T = 



pA^dx + i f 

11 2 J l t 


P i 'el dx + 1 m e :xf + 1 J ei ^ 2 




(15) 

where the first term describes that energy associated with the lateral motion of the cen- 
terline; and f refers to integration over the length of the i^ 1 beam. The second term, 
^i 

the rotary inertia term, defines the energy associated with the rotation of an infinitesi- 
mal slice with the angular velocity 6. The following two terms result from the transla- 
tion and rotation of the engine; X and /3 are the lateral and angular velocities, respec- 
tively. The last term is the contribution of the sloshing mass. 

Differentiating the displacement equations (eqs. (1), (2), (4), (6), (8), (9), (11), 

(12), and (14)) with respect to time and substituting the resulting equations into equa- 
tion (15) yields a quadratic expression in the generalized velocities. Expanding this ex- 
pression and interchanging the summation and integral signs provides the following 
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kinetic energy expression: 


f o f o 


= | ^ ^ ^^On^Op^+^^oVop^On'lOp 


p=T, R, 1 n=T, R, 1 

r f, 


+ I m ^ 


0 


E Kn( x eO> + e 0^0n( x e0^0n + e 0«0 | 

n=T, R, 1 


+ - J 


eO 


-.2 


E ^0n( x e0^0n + «0 


k / f 0 f O 


r u v ^ 

^ 2^ / y jjf pA i[^On^ + ( x - x i^OnMhp( x i) + < x - x i^OpM dX 


i=l \ p=T, R, 1 n=T, R, 1 


f i f O 


+ p VW x i^Op^ fc |%n p Op + V ^ (/ Mi Kn< x i> + < x " X i^OnM ^ip 

- < m ti -I ^ 


p=l n=T, R, 1 


f i f i 


+ f h Ph+oMn P ^Wn% + l mi P Vin<Pip ^ + l Pi^in^ip ^i^ip 


p=l n=l 


f i 


+ \ m eii E [^On( x i) + ( x ei + e i " x i^On< x i)] p On + E frin< x ei> + e i^i n < x ei>] ^in + e i«i 


I n=T, R, 1 


n=l 




f i 


-i 2' 


E ^On^On + E ^in( x ei^in + “i 


n=l 


v 

+ E \ “r E Kn( x r)%n + ^r ] ‘ 


n=T, R, 1 


( 16 ) 
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The kinetic energy expression can be simplified by introducing the following conditions: 
(1) Conservation of linear and angular momentum for the uncoupled elastic free-free 
modes of the main beam; that is, preservation of translational and rotational 
equilibrium 


J' PA-o^On ^ - ® (n - 1, 2, . . . , fg) 
l 0 

f P A o x <P 0 n + f P I 0^0n dx=0 (n= 1, 2, . . f Q ) 


(2) Inertial orthogonality and normalization condition for the uncoupled modes of the 
branch beams as shown in appendix A 


/ PV in^ip dx + / Ph+ in^ip ^ = 6 npPi 
L i 


(i — 0, 1, 2, . • . , k) 


(ft? P ~ 1? 2, • . . 9 f •) 


where 6 np is the Kroneeker delta defined such that 


0 (P ^ n) 


6np |l (P = n) 


and p in is the generalized mass of the n^ uncoupled normal mode of the i*"* 1 beam. 

Since the elastic modes of each branch are orthogonal relative to each other, a 
number of inertial coupling terms are eliminated from the kinetic energy expression. 
Those that remain are generally easily evaluated. To simplify the remaining terms 
further, the following notation is introduced: 


m 0 B 0 - J 

l 0 

pAgX dx 

(20a) 

= 1 PA 0 x 2 dx+ f pl 0 dx 
*0 *0 

(20b) 

pA^x - Xj)dx 

(i = 1, 2, . . . , k) 

(20c) 
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I 


I II Mill 


Jj = f pA t (x - x.) 2 dx + J' pl t dx (i = 1, 2, . . . , k) 




'*i 


Q. 

T.n 

~2 

co- 

in 


r /i - 1, 2, , k\ 

' •( PV ‘" ^ \n = lj 2, . . . , fj/ 


(20d) 


(20e) 


M in r r /i = 1, 2, 

= / PA^x - Xi V in dx + J pl^ in dx ( 

0)“ ^ i \ n_ L > *> 

in 


. k) 

• • , v 


(2 Of) 


where 

m 0 

B n 


“i 

B i 

J, 


Qi 


in 


O)- 

m 

M in 


mass of main beam 

distance from center of gravity of system to center of gravity of main beam 
mass moment of inertia of main beam about center of gravity of system 
mass of r branch 

distance from center of gravity of system to attachment point of i** 1 branch 
mass moment of inertia of i^ branch about attachment point 

modal shear at clamped end of i^ 1 branch in n^ 1 uncoupled cantilever mode (see 
fig. 3) 

natural frequency of n^ uncoupled cantilever mode of i** 1 branch 

modal moment at clamped end of i** 1 branch in n^ uncoupled cantilever mode (see 
fig. 3) 
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The simplified kinetic energy expression now appears in terms of total masses, mass 
moments of inertia, natural frequencies, generalized masses, and generalized forces 
(modal shear and moment at base) of the components: 


T “ o ( m 0^0T + 2 m 0 B 0^0T^0R + J cAoR + 2 M 0n^0n 


+ \ m e 01 X) [^On^eO* + ^On^eO^On + e 0 a 0 


n=T, R, 1 


+ 7 J e0 2 ^0n( x e0^0n + “0 


n=T, R, 1 


k / f 0 f 0 


2 U 2 2 K'POnW^Op^i) + m i B i^ 0 n( x i^ 0 p( x i) 


i=l \“ p=T, R, 1 n=T, R, 1 


^VonWopty + ^On^Op^Wtop 


f i f 0 


p=l n=T, R, 1 


r i f i 

Qi M - _ 4 ■ ^ o 

-9 OnW ~T + 't'OnW ~Y Wkp + 7 2 ^in^ii 

w ip w ipJ n=1 


l m ei ] 2 ton^i) + < x ei + e i “ x i^0nM 


n=T, R, 1 


2[^in( x ei) + £ i^in( x ei)]v + e i«i 


U 1 

\ J ei 2 Von^On + 2 *in( x ei>*in + «i 
z n=T, R, 1 n=l 


r=l “ n=T, R, 1 


+ 2 ! m r 2 < f , 0iS\y%n +2 ’: 


15 


Thus far, the theoretical development has been restricted to beams; the kinetic en- 
ergy and the orthogonality conditions are those of a beam. Suppose that the branch can- 
not be represented as a beam, but rather as a spacecraft composed of trusses and dis- 
crete masses. In this case, the final kinetic energy expression (eq. (21)) is still appro- 
priate provided that the displacement can be expressed as the superposition of the trans- 
lation and rotation of the branch as a whole, and of its cantilever normal modes. The 
mass properties and displacement functions in the kinetic energy expression (eq. (15)) 
then correspond to discrete points. The energy contribution involves coupling between 
the generalized coordinates describing the motion of the main beam and the branch. The 
orthogonality conditions again prevail to eliminate cross coupling terms between the gen- 
eralized coordinates of the branch. The remaining coupling terms are simplified by ex- 
pressing them in terms of the mass and mass moment of inertia of the branch, and the 
modal shears and moments developed at the base of the branch in the cantilever modes. 


Potential Energy 

The potential energy associated with this system of beams, engines, constraining 
springs, and sloshing masses is 

k 

"E 

i=0 



+ \ X) K Tj y o( x Tj) + \ ^2 K e m 6 O < x 0 m) + \ 2 K r^r ( 22 > 
j m r 

The first term describes that strain energy associated with beam bending, and the next 
term, the shear correction term, defines the energy associated with the shear slope 
$i - (9y^/3x). The ensuing term includes the energy stored in the engine spring. The 
following two terms define the energy stored in the translational and rotational con- 
straining springs, and the last term is the energy of the sloshing spring. 

Differentiating equations (1), (2), (4), and (6), where the prime denotes differentia- 
tion with respect to x, yields 
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0y O _ V'' , 

IT " Z, ^On q On 

m ti 1 


n=T, R, 1 


3x 


^ x 

XJ ^On^i^On + ^ > ^in q in 


n=T, R, 1 


n=l 


30 


l O 


3x 


~~ X/ ^On q On 
n=T, R, 1 

f i 

1 = S ^’inltn 


30 


3x 


n=l 


(23) 


and substituting equation (23) into equation (22) yields the potential energy expression: 


+ K 


e0“0 



dx 


+ f ( A s G >. 

^in " ^ in^in 

2 ^ 
dx + K ei af > 


1 1 

n=l 




-> 


r 

f 0 

l 2 


+ 2^ K T3 


n=T,R,l 


!0n 


+ 



f 0 

£ 

n =T,R, 


^On^O m^On 

1 



(24) 
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The potential energy may be simplified by introducing the following conditions: 

(1) The rigid body modes and their derivatives must not contribute to the bending 
and shear strain energy of the beams: 

^OT = ^OR = ^OT " = ^ OR ~ = 0 ( 25 ) 


(2) The stiffness orthogonality condition of the elastic uncoupled modes of the branch 
beams (appendix A) must be 



- «>ip)dx 


= 6 


2 

U. OK 

np p m m 


(i - 0, 1, 2, . . . , k) 

( n ? P = 1? 2, . . . , f^) 


(26) 


o 

where 6 is the Kroneeker Delta defined in equation (19), and ju. co. is referred to 
np m in 

as the generalized stiffness of the n 11 uncoupled mode of the i Ln beam. 

Again, the orthogonality condition, characteristic of the normal modes, proves in- 
valuable in eliminating cross coupling terms. Its merits are obvious in the following 
potential energy expression, which results from substituting equations (25) and (26) into 
equation (24). 


f i 


U = 




i=0 \n=l 


X ^On^Tj^On 
n=T, R, 1 


IS 


K, 


■0 m 


m 


X ^0tS x 9 m^ q 0n 
n=T, R, 1 


|X K ^r < 27 > 


Examination of the potential energy expression reveals that the terms pertaining to 
the branch beams, swiveling engines, and sloshing masses have reduced to the canonical 
form ^ a.q 2 of the quadratic, while those terms associated with the constraining 
1 1 1 / \2 

springs remain ( 2J a-q^ 1 . When this function is substituted into Lagrange’ s equation, 


the quadratic terms will produce cross products of generalized coordinates (static cou- 
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pling terms) in the equations of motion. Equation (27) reveals that static coupling is in- 
troduced only by the constraining springs, while the addition of terms for branch beams, 
swiveling engines, and sloshing masses does not result in static coupling. 


Equations of Motion 

The kinetic energy and potential energy equations ((21) and (27)) have been evaluated 
in terms of the generalized coordinates. It is now possible to introduce Lagrange's 
equation 


d_ 

dt 



+ iu = o 


aq t 


(28) 


The application of this equation supplies the M f Tt equations of motion, each of which cor- 
responds to a specific generalized coordinate. Because of the formidable appearance of 
the equations of motion, it is natural to bypass the present method of notation by the in- 
troduction of matrix notation. In matrix form, the Tt f M equations of motion appear as 
follows: 


[B] { q } + [A] { q } = {0} 


(29) 


where 

[B] inertia matrix describing coupling of various masses of system 
{ q } column matrix of generalized coordinates 

[A] stiffness matrix defining coupling effect of various stiffnesses of system 

Again, the complicated nature of the matrices, especially the inertia matrix, warrants 
special consideration in their presentation. Each matrix is decomposed into several 
elementary matrices, and the stiffness and inertia matrices are presented as sums of 
these matrices. The elementary matrices pertain to specific components, such as 
branches, engines, sloshing masses, and constraining springs. In this manner, various 
components may be added to or removed from the system at will, so that 

k 

[A] - £) (l A il + f A ei]) + J2 t A srl + ( A t 1 + [ A e ] < 30 > 

i=0 r 
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( 31 ) 


k 


[ B ]'E([ B i] + Pei]) + I>sr] 

i=0 r 


where 

[M, [BJ 


t A eiM B eJ ] 

[ A s r]>[ B SI J 

[Arp] , [Ag ] 


J.U 

stiffness and inertia matrices, respectively, for i cn branch 
(i = 0, 1, 2 . . . ) 

i.U 

stiffness and inertia matrices, respectively, for engine on i cn branch 

th 


stiffness and inertia matrices, respectively, for r sloshing mass 

stiffness matrices for translational and rotational constraining springs 

The formation of these matrices can best be shown by an example. Consider a system 
that consists of a main beam with an engine, two branch beams with engines, a sloshing 
mass, translational constraining springs, and rotational constraining springs. The dis- 
placement of the system is described in terms of the following generalized coordinates: 

rigid body coordinates 

three elastic coordinates of main beam 

rotation of engine on main beam 


q 0T’ q 0R 

q 01’ q 02’ q 03 
a n 


q H, q-^2’ ^13 three elastic coordinates of first branch beam 


a 


q 21’ q 22 


rotation of engine on first branch beam 

two elastic coordinates of second branch beam 

rotation of engine on second branch beam 

j displacement of the sloshing mass 

The position of the generalized coordinates in the column matrix is the same as the list 
just given; for example, the transpose of the column matrix appears as 


{ q } - [ q 0T q 0R q 01 q 02 q 03 a 0 q ll q 12 q 13 a l q 21 q 22 a 2 ^lJ 


The general procedure for the placement of the generalized coordinates in the column 
matrix is the same as stated previously; rigid body coordinates (qg-p qgp), main beam 
and engine coordinates fog-p q 0 2> • • 

( q ip qj2’ • • •’ q lf 1 ’ a V q 21> q 22’ 
and the sloshing coordinates ( Jfcp 


q 0f ’ a o)’ ^ ranc h beam and engine coordinates 

’ * » q kf k ’ » 

, -S^ r ). The elementary stiffness and mass 


q 2f 2 ’ a 2’ 


q kl’ q k2’ 


matrices are now presented. 
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Main beam. - For the main beam, the stiffness and mass matrices are 
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Engine on main beam. - With an engine on the main beam, the stiffness and mass 
matrices are 


[ B e 0 ] = 


q 0T q 0R q 01 q 02 q 03 "0 q ll 


[ A e 0 l 


eO 


m 


eO 


r 1 

[Transpose of column matrix] 

q 0T 

x e0 + e 0 


q 0R 

^Ol^eO^ + € 0^01^ x e0^ 


q 01 

^02^ x e0^ + e 0^02^ x e0^ 


q 02 

^03^ x e0^ + € 0^03^ x e0^ 


q 03 

1 € ° 


a 0 

0 


q ll 

0 


'*1 


+ J 


e0“ 


" 0 " 

[Transpose of column matrix] 

q 0T 

1 


q 0R 

*0l( x e0> 


q 01 

^02^ x e0^ 


q 02 

o 

<u 

^co 

o 


q 03 

-i 1 

>■ 

a 0 

0 


q ll 

. 0 , 


*1 
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Main beam. - For the main beam, the stiffness and mass matrices are 



Engine on main beam. - With an engine on the main beam, the stiffness and mass 
matrices are 


q 0T q 0R q 01 q 02 q 03 a 0 q ll * * • 


0 


0 


t A e 0 ] = 


0 


0 





r 1 

[Transpose of column matrix] 

q 0T 

x e0 + e 0 


q 0R 

‘Poi^eCp + e 0^01^ x e0^ 


q 01 

^02^ x e(P + e 0^02^ x e0^ 


q 02 

^03^ x e(P + e 0^03^ x e0^ 


q 03 

1 e 0 

b 

“0 

0 


q ll 

0 


*1 


r 0 

[Transpose of column matrix] 

q 0T 

1 


q 0R 

^0l( x e0> 


q 01 

^02^ x e0^ 


q 02 

^03^ x e0^ 


q 03 

4 1 

b 

a 0 

0 


q ll 

0 J 


*1 
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V 


First branch beam . - For the first branch beam, the stiffness and mass matrices 

are 


q 0T q 0R q 01 q 02 q 03 a 0 q ll q 12 q 13 a l 


' *1 


[A x ] = 


M ll w ll 


q 12 w 12 


m 13 w 13 


0 


[ B l] 


‘11 


12 


q 13 


+ < 


<?oi( x i) 

^02^ x l^ 

^03^ x l^ 


-Q n -Q 12 -Q 13 


2 2 
; 11 w 12 


4 


q 0T 

q 0R 

q 01 

*02 

q 03 




+ rrij 


+ 


+ m^B 


q 0T q 0R q 01 q 02 q 03 


^ 1 

x. 


^Ol( x l> 

<pQ2^ x l> 

^ 03 ^ x l) 

0 


[ 1 x x ^q^Xj) <p 02 (x 1 ) <p 03 (x 1 ) 0 


q 0T 

q 0R 

q 01 

q 02 

q 03 


0 ^ 
1 

^Ol( x l) 

* 02 < x l> 

* OS^l* 

0 


[ 1 X 1 ^02^1^ ^03^ x l^ 0 


° J 


q 0T 
q 0R 
q 01 
q 02 
q 03 


L i 


<Poi( x i) 

(pQ2^ x l> 

V 03^ x l^ 


[ 0 1 >^01^ x l^ ^02^ x l^ ^03^ x l^ 0 


r 0 
1 


V'oi( x i) 

^ 02 ( x l) 

^ ^ 03 ^ x l^ 
0 


[ 0 1 ^0l( x l) ^02^ x l^ ^03^ x l^ 0 


° J 


0 J 


q 0T 

q 0R 

q 01 

q 02 

q 03 


q 0T 

q 0R 

q 01 

q 02 

q 03 
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Engine on first branch beam. - For an engine on the branch beam, the stiffness and 
mass matrices are 


Pel] = 


q 0T q 0R q 01 q 02 q 03 a 0 q ll q 12 q 13 a l q 21 q 22 a 2 *1 


Pel] 


0 


0 


0 


0 


0 


0 


0 


0 


HfSj. 


0 


k el 


0 


0 


0 


0 


r 1 

1_ Transpose of column matrixj 

q 0T 

x el + € l 


q 0R 

<?0l( x l> + ( x el + e l " x 1^0l( x l) 


q 01 

<P02^ x l) + ( x el + € 1 " x 1^02^ x l^ 


q 02 

^03( x l) + ( x el +e l ' x l>^03^ x l> 


q 03 

0 


“0 

*ll< x el> +e l*ll< x el> 

y. 

q ll 

^l 2 ( x el) +€ 1^12( x el) 


q 12 

^13^ x el^ + € 1^ 13^ x el^ 


q 13 

e l 


"l 

o 


-S' 

V- J 


1 


' 0 " 

[Transpose of column matrixj 

q 0T 

1 


q 0R 

^01 (X 1> 


q 01 

^ 0 2( x l) 


q 02 

*03< x l) 


q 03 

0 


"o 

J ^ll( x el> 


q ll 

^l 2 ( x el> 


q 12 

*13< x el> 


q 13 

1 


“1 

* 

. 0 ^ 
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q 0T q 0R q 01 q 02 q 03 «0 


+ m 2 « 


•Poi^ 
^02^2) 
*P 03^ x 2^ 
0 


[ 1 x 2 <PQ\(*2> *P 02^ x 2^ ^03^ x 2^ 0 


°J 


q 0T 
q 0R i 
q 01 
q 02 
q 03 


a 


0 


+ m 2 B 2 < 


^ 01 ^ x 2 ^ 

^ 02 ^ x 2 ^ 

^ 03 ^ x 2 ^ 

0 


l 1 x 2 ‘Poi^ ^ 02 ^ ^ 03 ^ 0 


°J 


q 0T I 
q 0R | 
q 01 
q 02 
q 03 
0 


a 




+ m 2 B 2 ^ 


x 2 

^ 0 l( x 2 ^ 

^ 02 ^ x 2 ^ 

03 ^ x 2 ^ 


[ 0 1 ^02^ ^03^ 0 


°J 


q 0T 

q 0R 

q 01 

q 02 

q 03 

«0 




+ Jo 


^ 02 ^ x 2 ^ 
J ^03^ 


[ 0 1 ^01^ x 2^ ^ 02 ( x 2^ ^ 03^ x 2^ 0 


q 0T 

q 0R 

q 01 

q 02 

q 03 

“0 


jgr 
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Engine on second branch beam. - For an engine on the branch beam, the stiffness 
and mass matrices are 


q 0T q 0R q 01 q 02 q 03 a 0 q ll q 12 q 13 01 1 q 21 q 22 a 2 


[ A e2^ = 


e2 


[ B e2l = 


= m 


e2 


r 1 1 

[Transpose of column matrixj 

q 0T 

x e2 + e 2 


q 0R 

^01^2^ + ^ x e2 + e 2 _ *2^01^ 


q 01 

^02^ x 2^ + ( x e2 + e 2 _ x 2^02^ x 2^ 


q 02 

<P 03^ x 2^ + ^ x e2 + e 2 " x 2^03^ x 2^ 


q 03 

0 


"0 

< 

y 

• 

<?? 21^ x e2^ + e 2^21^ x e2^ 


q 21 

^22^ x e2^ + e 2^22^ x e2^ 


q 22 

e 2 


a 2 

0 

t. J 


1^1 


+ J 


e2 


'' 0 ^ 

[Transpose of column matrixj 

q 0T 

1 


q 0R 

^0l( x 2^ 


q 01 

^ 02^ x 2^ 


q 02 

^03^ x 2^ 


q 03 

0 


a 0 

■< 

> 

• 

^21^ x e2^ 


q 21 

^22^ x e2^ 


q 22 

1 


a 2 

^ 0 j 


*1 
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